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Abstract 

We study a chain of coupled nanomagnets in a classical approximation. We show that the 
infinitely long chain of coupled nanomagnets can be equivalently mapped onto an effective one- 
dimensional Hamiltonian with a fictitious time-dependent perturbation. We establish a connection 
between the dynamical characteristics of the classical system and spin correlation time. The decay 
rate for the spin correlation functions turns out to depend logarithmically on the maximal Lyapunov 
exponent. Furthermore, we discuss the non-trivial role of the exchange anisotropy within the chain. 
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I. INTRODUCTION 



Nanoscale magnetic structures have promising applications as basic elements in future 
nanoelectronics devices and are frequently discussed in the context of quantum information 
processing. The principal challenge of the quantum information technology is finding an 
efficient procedure for the generation and manipulation of the many-qubit entangled states. 
Those can be realized on the basis of e. g. Rydberg atoms located in optical quantum cavities 
2-4], Josephson j unctions j^J, or ion traps 6]. One very promising realization is based on 
single molecular nanomagnets (SMMs)[7H9|. These are molecular structures with a large 
effective spin. A prototypical representative of this family of compounds is, for example, 
Mn 12 acetate in which S = 16. Molecular nanomagnets show a number of interesting 
phenomena that have been in the focus of theoretical and experimental research during the 



last two decades. 
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18] For instance, SMMs show a bistable behavior as a result of the strong 



uniaxial anisotropy |7| as well as a tunneling of the magnetization. |8| An attractive feature 
for information storage is the large relaxation time of molecular nanomagnets. 12J 

SMMs are usually modelled by spin chain Hamiltonians augmented by different kinds of 
interaction terms responsible for different compounds. These interaction contributions are 
highly non-trivial and in most cases are anisotropic. This makes the analytical treatment 
very cumbersome calling for efficient theoretical approaches. In this paper we investigate 
the properties of a classical spin chain coupled by anisotropic exchange interactions. This 
case is relevant not only for chains of exchange coupled SMMs, 2l|] but also for several other 
realistic physical problems. Those include weakly coupled antiferromagnetic rings 
large spin multiples coupled by Dzyaloshinskii-Moriya (DM) exchange interaction 
will demonstrate that even for multidimensional complex physical systems it is still possible 
to obtain analytical results using special mathematical technique presented in 24J. Its 



,8] 




applicability to the SMMs is shown in 
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231 ] . In the first step one evaluates the Lyapunov 



exponents and the spin correlation functions for the system. Having done that one can then 
extract information on the properties of the system beyond the classical limit. For example, 
there is a deep connection between the classical Lyapunov exponents and the quantum 
Loschmidt echo 25|, which is a natural measure of the quantum stability and of the fidelity 
of quantum teleportation l|. If the Lyapunov regime is reached for a quantum system, then 
the decay rate for the teleportation fidelity can be identified via the Lyapunov exponent. 
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Formal criteria for the Lyapunov regime 26| in case of a chain of SMM can be estimated from 
the relation A < J 2 / A, where A is the Lyapunov exponent, A is the mean level spacing and 
J is the exchange interaction constant between SMM spins (which prohibits the integrability 
of the system leading thus to the chaotic dynamics). Furthermore, it can be shown that the 
spin correlation functions can be expressed through the Lyapunov exponent. 

The paper is organized as follows. In Section II we give a brief exposition of our model 
and present the details of the principal investigation technique. After that in Section III 
we discuss the relation between the relevant spin correlation functions and their classical 
analogs. Section IV contains the treatment of the system in the chaotic domain. Finally, 
the conclusions section summarizes our findings. All necessary mathematical details are 
presented in Appendices. 

II. THEORETICAL FORMULATION 

The prototype model Hamiltonian for the exchange coupled SMM is 

h = j SnS: l+1 + g £ (s*s* n+1 + sy n si +1 \ ( s n) » (!) 

where J and g are exchange interaction constants, = —DS 2 stands for the anisotropy 
barrier height of the system. For the prototypical Mriyi acetates |tJ D ~ 0.7K sets the 
value of the barrier parameter (23] and S x ' y,z are spin projection operators of SMM. Due to 
the large spin of the SMM (see Introduction) analytical quantum mechanical treatment of the 
model (CQ) is hardly accessible. To make progress we choose the semi-classical parametrization 
as follows: 

S z n = cos 9 n , SI = sin 9 n cos <p n , = sin 9 n sin yv (2) 
Then the Hamiltonian (P) can be rewritten in the more convenient form 

H = J ^2 cos 9n cos 9 n+ i + g £ sin 9 n sin 9 n+1 cos(</?„+i - <p n ) + f3 £ cos 2 9 n . (3) 

n n n 

Our aim is the evaluation of the correlation functions and the study of the spin dynamics 
governed by the Hamiltonian ([3]). Since this is highly nonlinear problem, it cannot be done 
in a simple and direct way. However, one can rigorously show, that there is a direct map 
between the chain of SMM and a one- dimensional (ID) model Hamiltonian with a fictitious 
time-dependent external perturbation. 



The equilibrium state for the model (3) satisfies the minimum condition of the infinite- 
dimensional functional H[9,(p]: 

dH dH 

o9 n d(f n 
n = 1,2... oo. 

Considering the Hamiltonian ([3]) we retain only the first order terms of the anisotropy 
parameter e = ( J — g)/2g <C 1. Then, after straightforward but laborious calculations we 
deduce from (Tj0) that the following relations hold 

S n+1 = (-l) m [S n - /3sin(20 n )(l - ecos(20 n ))], (5) 
9 n+l = (-l)™e n + 7rv+ (-1)" arcsin S n+1 , (6) 
<p n+ i = ip n + Trm, m = 0,1 v = 0, 1, 

where 5* n+ i = sin(# n+1 — 6 n ). The index z/ refers to the two possible solutions when inverting 
the trigonometric expression for S n+ i. Depending on the sign of the re-scaled barrier height 
— > =F (3/g the index m = 0, 1, defines the energy minimum condition. For convenience we 
will use positively defined > and consequently m = in our calculations. The above 
result is just a recurrence relation in the form of the explicit map [S n+ i, 9 n +i) = T(S n , 6* n ). 
Our idea is to find a Hamiltonian model which is equivalent to (JSJ). Let us consider the 
following perturbed ID Hamiltonian system: 



H = H (s) + f3V(6)T S(t-nT), 

n=—oo 

H (S) = vnS + (-I)" (s arcsin S + VT^P 



(7) 



V{9) = -^cos 2 -|cos 2 29 

The respective Hamiltonian equations read 

dS OH 



= -PV\9)T 6(t-nT), 

n=—oo 

u(S), (8) 



dt 89 

d9 _ dH 
~dl~~dS " 
ui v {S) = itv + (—I)" arcsin s 



Their integration simplifies due to the presence of the delta function in the perturbation 
term because the evolution operator (S; 9) = T(S; 9) splits in a pulse-induced Ts and free 



evolution terms T R : 

f = f R ® f 5 , 

S = S(t Q + t-0); 9 = 9{t + T -0), (9) 

s = s(t -oy, e = e(t -o). 

For operator of free evolution we simply have 

f R (S;9) = (S;9 + u u (S)T). (10) 

The explicit form of the pulse-induced evolution operator T$ can be derived after an inte- 
gration of the system (6) on the small time interval (to — 0, to + 0) around t where the pulse 
is applied 

*o+0 t +0 

S(t Q + 0) - S(t - 0) = 



J Sdt = - J P^p-T 5(t-kT) = -(3T 



dV{9) 



io-0 t -0 
to+0 



89 

k=—oo 



0(t o + 0) - 0(t o - 0) = J 9dt = 0. (11) 

to-0 

Taking into account Eq. (fTTI) . for the pulse- induced evolution operator we obtain 

M s -^4 (12) 

Combining Eqs. ( TTUj) and (fT2|) the complete evolution picture can be expressed through the 
following map 

(S, 9) = f(S, 9) = f R f s (S, 9) = f R (s - PT 9V< & 6 \ ej = 

S-P?^p-;9 + u v (S)y (13) 
Or in the explicit form: 

S n+1 = S n - (3 ( sin 29 n - | sin A9^j , 

n +i = n + u v (S n+1 ), (14) 
Uv(S n ) =txv + (-l^arcsinSVi, 

For the details of the derivation of Eq.(14) see Appendix A. 



This result is obtained for the kicked Hamiltonian model with T = 1 and it matches 
exactly the recurrence relations in Eqs. (jSJ) obtained for the SMM chain. Such an analogy is 
quite important since the infinite-dimensional nonlinear system is now equivalent to the 
ID Hamiltonian model. We note, the discrete time in the perturbation term (j7j) is fictitious 
and corresponds to the number of the spins in the chain (pQ). 



III. SPIN DYNAMICS AND CORRELATION FUNCTIONS 



We will proceed with the equivalent ID Hamiltonian model with fictitious time dependent 
perturbation (jTJ), which is more convenient than the multidimensional nonlinear model (j3]). 
Due to the nonlinearity of the model (ED) we should expect a rich and a complex dynamics. 
In particular, our purpose is to establish a connection between the chaotic dynamics and 
the decay rates of the spin correlation functions. As a first step we construct the Jacobian 
matrix of the map (TfJ 



(as as \ 




as ae 


H 


1 89 88 
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1 pV{d) 
u'(S) 1- (3u'(S)V"(9) 
The Lyapunov exponents can be evaluated as eigenvalues of this matrix and are thus given 
by 

Ai,2 = i + f ±^(i + ljr) a -i, (16) 

where 

K = -/3u;' v (S)V"(9) (17) 
is the chaos parameter 271 ] . From Eqs. ffTBl) and (|17[) for the chaos parameter K we find the 



following simple relation 

2/3 (cos 20 -e cos 40) , . 

K = — . '-. 18 

The dynamics is expected to be chaotic if K > 0, Ai > 1 or K < —4, A2 < — 1. Therefore 
from Eq. f fl8|) we obtain the relevant intervals for the angle variable 



9 E ( ivm H — arccos ( — 
V 2 \2e 

m = 0,±l.... (19) 



1 k* +2 )> {m+1) *-h TCC0S (h-y^ 



In the isotropic case J ~ g, e — > 0, this leads to 



9 E [ nm + j; (m + 1)tt - - ), m = 0, ±1 . . . . (20) 



-1.0 



FIG. 1: Color online. The parameter K(9,S) signifying chaotic behavior plotted according to Eq. 
(|T8|) for the following values of the parameters [3 = 3,e = 0.5. 

Equation (19) defiance the width of the chaotic domain where parameter of chaos (18) is 
larger than one K > 1. Obviously, the width of the chaotic domain depends on the values of 
anisotropy and easy to see that in case of small anisotropy area of chaos is narrower than in 
zero anisotropy case A8(e = 0) > A6(s << 1). Therefore we conclude that small anisotropy 
leads to the less chaotic regime. 

From the parameter of the chaos (IT8"|) (also plotted in Fig. [T|) we conclude that the phase 
space of the system consists of domains corresponding to a regular and a chaotic motion. 
Later we will use the random phase approximation, which is valid precisely in the latter 
domain, j^] 

In order to obtain explicit expressions for the spin correlation functions we rewrite the 
recurrence relations in Eqs. ( |T4l) in the following form 







+1 = 6 n + u(s n+1 ) = 9 n + uu u (S n - pv'(e)) 

e n + u v {s n )-i3v\e)uj l u {s n ) = 



(21) 




For the angular variable we infer the self-consistent recurrence relation 




(22) 
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The correlation function is given by 



2- 



Ztx J 
o 



2tt 



1 

2^ 



'n-l+avfSn-l)— -Ko sin26» n _i + — ^ sin40„_i — 0o 



(23) 



and can be calculated using the above iterative procedure as well as the expression for the 
Bessel function exp(iz sin ip) = J m (z)e mv . Taking into account that u v (S n ) = u ~ 

const, Kq = 0cu'(S n ) ~ const, from Eq. f[2"3"j) we deduce 

2tt 



2tt 



1 



y^Q J{9„-i-9o) -iK sin28 n - 1 ^ £ sin46 n _i 




271 J 







il=— oo h=— oo 



(24) 



hoo +oo 



-oo +oo +oo 



+oo 



+l 2 + ...+ln 



= «*" E E ••• E E E ••• E 

mi=- oo m,2=— oo m n =— oo /i=— oo l2=— co £ n =— oo 
2iuJraig— 2iu;(mi+m2) g-22cj(mi+m 2 +...+7n n _i) e -2io;/i e — 2iw(Zi4-/2) g— 2iw(Zi+Z2+— +2n-i) 

•J mi [A^ ]«/ m2 [(l - 2mi + 4h)K ] . . . J m „[(l - 2(mi + m 2 . . . + m n _i) + 4{h + l 2 ... + l n ^))K Q } ■ 



K e 



. . . Jl 



1 - 2(mi + m 2 . . . + rvi) + 4(Zi + Z 2 • • • + l n -i))K e 



•01-2(mi+m2+...+m„)+4(Zi-H2+---+Z«);l ■ 

For the details of the derivation of Eq.(24) see Appendix C. 

In the case of a large Lyapunov exponent J m (K ) ~ (Kq)' 1 / 2 , K ^> 1 we infer from ( 124 



exp 



n 



(25) 



where r c = 2/ In 
estimation 



, is the correlation length. Since K ~ /3, we have the following 

2 



T, 



ln(^ 



(26) 
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In the isotropic case e — 0, one can perform the same calculations [insertion of e = into 
Eq. fl26|) gives a wrong result] and show that 



holds. Taking into account Eqs. (I26I) . ( )27|) and expressions for the rescaled interaction 
constants e — > B - 

Namely, the strong anisotropy 



constants e — >■ /3 — >• ^ we conclude that the role of the anisotropy is not trivial 



J -9> A -J, (28) 
suppresses the spin correlations because then r' c > r c . On the other hand the weak anisotropy 

J -9<^, (29) 

enhances the correlations t' c < r c . 

It should be stressed that the reliability of the analytical estimates is limited. This is 
particularly apparent for the case where the numerical and the analytical predictions deviate 
from each others, due to limited range of applicability of the analytical expressions, derived 
after rough approximations. 

The role of the anisotropy e — (J — g)/2g can be clarified numerically as well. In order 
to better understand the physical features of the model (CEJ) we will study the phase portrait 
of the system. The results of the numerical evaluation of the recurrence relations ( IT4|) are 
presented in Figs. [2j|5j As we see from Fig. [2] the phase space of the system consists of two 
topologically different domains separated from each other by a separatrix. Most of the phase 
space belongs to the domain of the regular motion and open phase trajectories. The domain 
of closed phase trajectories mainly corresponds to the irregular motion and a small island of a 
regular motion is observed only in the center of the portrait. From this formal mathematical 
statement, one can extract interesting physical information. Closed trajectories belong to 
the oscillatory regime and open trajectories to the rotational one. Therefore, we expect 
that two types of motion can be realized for the model ([1]) . The domain of the regular spin 
rotational motion is defined by the relation 0.2 < \S n \ < 1. If then \S n \ < 0.2 the spin 
oscillation is chaotic and only for very small amplitude an island of regular oscillations is 
observed in the center. If the anisotropy parameter is zero e = then the island of the regular 
oscillatory motion disappears (see Fig. [3]). This means that without the small anisotropy the 
spin system is less correlated [see Eqs. f )28|) and f l29|) ]. Such a geometrical interpretation can 
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Tt/2 71 
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FIG. 2: Color online. Results of the numerical calculations of the recurrence relations (|14p on the 
phase plane (S n ,9 n ). The following the parameters are used f3 = 0.05, e = 0.005. About hundred 
trajectories are generated for the set of different initial conditions (Sq,0o). 

be extrapolated from the pair of the canonical variables (S n , 9 n ) to the real spin variables 
S* = cos9 n using the parametrization (|2J) and a simple relation S n = sm(9 n — re _i). 



IV. SPIN DIFFUSION AND KINETIC APPROACH 

The dynamical picture does not apply in the chaotic regime for K > or K < —4. 
An adequate language in this case is the statistical approach. Instead of the dynamical 
variables the key role is played by the probability distribution function, which is a solution 
of the Fokker-Planck equation. Its derivation is rather straightforward for chaotic dynamical 



models and is based on the Kolmogorov, Arnold, Moser (KAM) theory. 24J Interested reader 



can find all technical details of the derivation for the spin chain model in the recent work [25]. 
Here we are using the final result adopted to the SMM system. The probability distribution 
of the spin variable S n = sin(6 l n — n _i) is described by following diffusion equation: 

df d 2 f 
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FIG. 3: Color online. Results of the numerical calculations of the recurrence relations (|14p on the 
phase plane (S n ,9 n ), for following the parameters: j3 = 0.05, e = 0. About hundred trajectories 
are generated for the set of different initial conditions (So, #o)- 

where D(S) = ^ ^1 + is the diffusion coefficient. For the details of derivation of Eq.(30) 
see Appendix B. The fundamental solution of this equation is 

'^db-K-ss)- (31) 

and can be found in many classical textbooks (see e. g. 30] ). This solution f[3"Tj) is defined on 
the interval — oo < S < oo whereas we need one for — 1 < S < 1. In order to find a solution 
relevant to our problem we will consider the following boundary and initial conditions for 
the diffusion equation (130]) . 

f = W for t = 0, 

f = gi{t) for S = -l, (32) 
/ = g 2 (t) for S = l, 

and will look for the solution in the following form 

oo 

f(S,t) = 2 sin(mTrS) exp ( - Dm 2 TT 2 t) M m (t), (33) 

m=l 
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ji/2 



3ti/2 



FIG. 4: Color online. Results of the numerical calculations of the recurrence relations (114D on the 
phase plane (S n , n ). We used f3 = 0.3, e = 0.15. About hundred trajectories are generated for the 
set of different initial conditions (So,9q). 



where 



M m (t) = J / (0 M™£)<% + D ™ J ex P {DmVr) [ 9l (r) - (-l) m g 2 (r)]dr. (34) 



In the simple case gi(f) = g2(t) = we obtain from Eqs. (I30l)-(l34 
AW 



f(S,t) 



E fo \u sin (( 2m + X > S ) ex P ( " D{2m + 1 > H ) ■ 

' 2m + 1 7 v 7 



7T ^ (2 

m=0 v 



(35) 



+1 +oo 

where the coefficient Wq can be denned from the normalization condition J J f(S,t)dSdt 

-1 o 

1, W = YTs)- Here <?(s) = Yl k~ s * s the Riemann zeta function 
^ ' k=i 



2s 



We note the direct 



correspondence between the fictitious time and the spin index t — > nT, T = 1 Eq. 
For the averages of the discrete random variable S = S n = sin(# n — 6 l n „ 1 ) we follow the 
standard procedure (see [22J, Eq.(18)-(20)) and utilize the distribution function Eq. fl35|) . 
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FIG. 5: Color online. Results of the numerical calculations of the recurrence relations fll4f) on the 
phase plane (S n ,6 n ), for the following parameters: /3 = 0.3, e = 0. About hundred trajectories 
are generated for the set of different initial conditions (Sq,0q). With the increase of the value of 
parameter (3, the domain of the chaotic motion covers almost the whole phase space (See Fig. [U 
Figj5]), because the chaos parameter K Eq.(18) is proportional to the constant f3. 

The integration is performed over the interval — 1 ^ S n ^ 1. As a result we obtain: 

+i 



< S 2 >=< S 2 >=< sin 2 



AW, 



^-i)»= J S 2 f(S,t)dS 



_o 

7T 



If 

V- / S 2 sm((2m + l)TT\S\)dS -exp(- D(2m+l)n 2 t). 

*-^ n (2m + 1) J 

m=0 v — 1 



(36) 



After an integration we get: 

< S 2 >=< sin 2 (0„ - dn-x) > 

7T 4 ~ V '~'2 



^e- DA $(e- 2DA ,4,i), (37) 



oc 



k=0 



where F({a\ . . . a p }; {b\ . . . b q }; z) = Yl WhTW)k ' ^ * S ^ e S enerai i ze d hypergeomet- 



ric function and $(z,s,a) = ( fc + a ) S is the generalized Riemann zeta function. [28 1 

k=0 

From Eq. (|37|) we immediately see that the statistically averaged random alignment fac- 
tor <C sin 2 (6> n — 6 n _i) ^> is not uniform along the spin chain (see Fig. |6|), but rather decays 
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FIG. 6: Results of the numerical integration for the statistically averaged random alignment factor 
■C sin 2 (# n — 9 n -i) S> in the diffusive approximation (|37p. for /3 = 0.3 and e = 0.15. 

exponentially with n. This result is reasonable since the solution for the distribution func- 
tion (135]) is obtained via deterministic initial and boundary conditions fl32|) . Therefore, a 
maximum of correlation is expected for n = 0. Since t = nT, t = 0, n = corresponds to the 
boundary where the distribution function is defined precisely. Far away from the boundary 
that means for n » 1 randomness occurs and the correlation decays. 



V. CONCLUSIONS 



In this paper we considered an anisotropic nonlinear spin chain, which serves as a model 
for a chain of coupled nanomagnets. We have shown, that there is a direct map between 
an infinite-dimensional spin chain model and an equivalent effective ID classical Hamilto- 
nian with a discrete fictitious time-dependent perturbation. We have established a direct 
connection between the dynamical characteristics of the classical system and the spin corre- 
lation time of the original quantum chain. The decay rate for the spin correlation functions 
turns out to depend logarithmically on the maximal Lyapunov exponent. In addition, for an 
anisotropic couplings we found an interesting counterintuitive feature: the small anisotropy 
leads to the formation of small islands of the regular motion in a chaotic sea of the system's 
phase space. As a result, the spin correlations become stronger within the islands of regular 
motion. We argue that these results obtained within the classical approximation are inter- 
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esting in other regimes. If the Lyapunov regime is reached for a quantum system, which 
takes place for the Lyapunov exponent A < J 2 / A, where A is the mean level spacing and J 
is the exchange interaction constant between spins, then the decay rate for the teleportation 
fidelity in a device based on such spin chains is directly related to A. 
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the Grant No. KO-2235/3, and STCU Grant No. 5053 is gratefully acknowledged. 

Appendix A: Derivation of the recurrence relations 

Let us consider the equilibrium state for the model (J3]): 

dH ( \ 

cos(6> n - 9 n+x ) - cos(6> n + 9 n+1 ) ) sin(y2 n+ i - tp. 



difn 

- 1 cos(6 l n _i - 9 n ) - cos(6» n _i + 9 n ) ) sin(y? n+ i - <p n ), (Al) 



dH J 



d6 n 2 



sin(0 n + 9 n+1 ) + sin(0„_i + 9 n ) + sin(0 n - 9 n+x ) - sin(0 n _i + 9 n ) - 



~ ( sm{9 n - 9 n+1 ) - sin(6» n _i - 9 n ) - sm{9 n + 9 n+1 ) - sin(8 n - 1 + 9 n ) ) cos(^ n+ i - ip r 



2 

+(3sm29 n = 0. (A2) 
After introduction of the notation S n = sm(9 n — 9 n _i), from flAip . ( 1A2j) we find: 

M) 5 --(M) s »-(M) si ^ +e »«)- 

I _ | \ S in(^ n _ 1 + 9 n ) + (3 sin(2^) = 0, (A3) 
(p n+ i = (f n + vrm, m = 0, 1 . 
Let the asymmetry parameter be defined by e = \J — g\, e < J,g. Next we perform a 



rescaling of the interaction constants ^ — > e, f3 — > =p^-. From flA3j) we deduce: 

(5^1 -S n )-e sin(9 n + 9 n+1 )e - sm{9 n ^ + 9 n ) + f3 sin(2^ n ) = 0, (A4) 
9 n+1 = 9 n + ixv + (-1)" arcsin^+i], v = 0, 1 . 

Depending on the sign of the rescaled barrier height /3 — > =p^, the value of the index m = 0, 1 
defines the energy minimum condition. For convenience we will use positively defined /3 > 
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and consequently m = 0. In the simplest case v = 0, so that from (1A4j) we obtain: 

(5*71+1 — S n ) 



-e^n[26 n ]y/l - S 2 n+1 + cos(29 n )S n+1 + sin(20 n ) ^1 -S*- cos(29 n )S n+1 j + (A5) 
+P sin(2# n ) = . 

Retaining only the first order terms with respect to the small parameter e = 2g , from 
( 1A5I) we obtain the following recurrence relations ( |T4|) : 

S n+1 = S n - (3 ( sin 29 n - | sin 46^ n ^ , 

#n+l = #n + ^(S'n+i), (A6) 

cu v (S n ) = ixv + (— l) v arcsin S n . 
Appendix B: Derivation of the kinetic equation 

The starting point for the derivation of the kinetic equation is the equivalent effective 
Hamiltonian ((7j): 

dS dH 



= -PV'{e)T Kt-nT), 

71= — OO 

u(S), (Bl) 



dt d6 

d6 _ dH 
~dl~~dS " 
u) u (S) = ixv + (— l) v arcsin s . 



Here the variable S plays the role of the adiabatic (slowly varying) action variable, while the 
angular variable 9 is the fast variable. Due to the presence of the two different time scales 
in the system: 

H = H (S) + eV(S,6,t) (B2) 
for the derivation of the kinetic equation we will follow to the standard procedure. 29 1 



The distribution function of the random variable f(S,t) obeys the Liouvillian equation of 
motion: 

U = iu(S)^ (B3) 



Li = -i| 



dV d dV d 
~dS86 ~ ~dSd6 j 
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The formal solution of the Liouville equation with the accuracy of second order terms in the 
small parameter e reads: 

t ti 

/DOS',*) = f(S 7 0) -ieYsf dtiexp[im j uj{t')dt']{n\L x \m)f Q {S,V>) + 



o 



+ (-ie) 2 J2 J dt i J dt 2 exp[-im J w(t')^ / ](0|L 1 (t 1 )|m)(m|L 1 (t 2 )|0)Zo(^ 0). (B4) 



Here (n\Li\m)^ J d9e~ mt 1 L\t %md is the matrix element of the Liouville operator. After 

n o 

averaging over the initial phases f(I,t) =<^C fo(I,t) 3> and applying the random phase 

approximation with respect to the fast chaotic variable ^(t^^i) = J u(t')dt' = 9{t{) — Qfa) 

t\ 

< expim\P(t 2 ,*i) >~ exp(— (ti - t 2 )/r c ) exp(-imoj(ti - t 2 )) ■ (B5) 
From (IB4p we obtain: 

9f_ = _ 2 2 v v (V^KQl^pMHL^pio^ fR s 

^ Z^Z^ (l/r c ) 2 + (mu-pO) 2 ' 1 J 

m>0 p>0 v ' y v y ' 

where 

(0|L lp |m)(m|L 1 _ p |0) = 

T/2 T/2 2tt 

^)lhf ] dtl J ^/^/^^(^^-^(^e-* ^ (B7) 

-T/2 -T/2 

and the following notation is used r c = 2T/lnK, T = After calculating the integrals in 
(1B7[) . in the limit ^ > 0, T = 1 form (1B6I) we simply recover the diffusion equation (1301) : 



W1 -^W, (B8) 



dt OS v ' <9S 

#4 



More details of derivations can be found in Ref. 
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Appendix C: Correlation functions 



For the evaluation of the multiple series in Eq. 
contributions from the main non-oscillatory terms. 

17 



(124"]) one should sum up the 
Due to the delta function 



Si- 2 (m 1 +m 2 +...+m n )+i(i 1 +i 2 +...+in),i in E( h (ED, and the fast exponential factors e iMh+h+-+i n ) ^ 
e -2io;(mi+ma+...+mn) ) ^he re l eva nt terms in Eq. fl2D are those with 

mi + m 2 + . . . + m n = 0, 

k + l 2 + ... + l n = 0. (CI) 

Using the asymptotic expressions for Bessel functions: 

J m (K ) ~ K " 1/2 , for A > 1, (C2) 
and condition (IClj) . one can easily obtain (26) from (T2] 
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